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We have computed the Bethe logarithms for the 1 1 S, 2 1 S and 2 3 S states of the helium atom to 
about seven figure-accuracy using a generalization of a method first developed by Charles Schwartz. 
We have also calculated the Bethe logarithms for the helium-like ions Li+, Be++, O b+ and S 14+ 
pH , for all three states to study the 1/Z behavior of the results. The Bethe logarithm of H - was 

O j' also calculated with somewhat less accuracy. The use of our Bethe logarithms for the excited 

states of neutral helium, instead of those from Goldman and Drake's first-order 1/Z-expansion, 
reduces by several orders of magnitude the discrepancies between the theoretically calculated and 
experimentally measured ionization potentials of these states. 
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PACS numbers: 31.15.Ar, 31.30.Jv 



Ever since the invention of quantum mechanics, the helium atom has served as an important testing-ground for 
our understanding of fundamental physics. In 1929 Hylleraas' calculation of the binding energy of the non-relativistic 
helium atom Hamiltonian showed that Schroedinger's formulation of quantum mechanics provided a quantitatively 
accurate description of not just two-body but three-body systems [jlj. During the 1950's, with the advent of fast 
digital computers, calculations by Kinoshita M and Pekeris || of not only the non-relativistic binding energy but 
also of relativistic corrections of 0(cv 2 ) Rydberg greatly improved the agreement between theory and experiment, 
and showed that the estimation of 0(a 3 ) Rydberg effects arising from quantum electrodynamics was important for 
obtaining agreement between theory and experiment at the level of 1 part in 10 6 or better. During the 1960's and 
1970's the variational techniques employed by Pekeris on the lowest states of singlet and triplet symmetry were 
extended to a wide range of excited states of the helium atom ||] During the 1980's, with the advent of two- 
{SJ , photon spectroscopy with counterpropagating laser beams, which can be used to eliminate the lst-order Doppler 
shift due to the thermal motion of the atoms, it became possible to measure the wavelengths for transitions between 
excited states of the helium atom with a precision of 1 part in 10 9 or better ||. Though numerous examples of 
excellent agreement between theory and experiment in a wide variety of contexts leave no reasonable doubt that 
C/2 | quantum electrodynamics is the correct theory for describing the interactions of charged particles at low energies, 
the extraordinary accuracy recently achieved in high-precision measurements on the helium atom poses a challenge 
to theorists to develop computational techniques capable of matching such accuracies. Since a 3 is of order 10~ 6 , it is 
clear that the coefficient of the lowest-order QED corrections needs to be evaluated with a relative accuracy of 10 
or better, and the effects of contributions with higher powers of a must also be estimated, to match the experimental 
accuracy of 1 part in 10 9 or better. 

For a helium atom or helium-like ion of atomic number Z, the leading 0(ct 3 ) Rydberg contribution to the Lamb 
| shift is given by the expression Q| 
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where the so-called Bethe logarithm H is defined by an infinite and slowly-convergent sum over all bound and 
continuum eigenstates: 

ln(fc /Ry) = | = 
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Here p is the sum of single-particle momentum operators (p = J2iPi) an d is an eigenfunction with eigenvalue Eq 
of the Hamiltonian H of the atom. For simplicity, we assume that H is the nonrelativistic Hamiltonian of an atom 
with atomic number Z, with a point nucleus of infinite mass: 



H = T+v = j2pin-zJ2 l /^+J2 1 / 7 
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which has the important and useful property that it is unitarily equivalent to the scaled Hamiltonian 
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which after division by Z 2 tends to a well-defined limit as Z — > oo. (The effects of the reduced mass /i = m e Mjv / (m e + 
Afjv) due to the finiteness of the nuclear mass Mn are subsequently included by scaling by appropriate powers of n/m e , 
and the negligible effect of the 'mass-polarisation' term M^ 1 Y^i>j Pi 'Pi on the Bethe logarithm is here ignored.) With 
the help of the closure relation J2 n \ n ){E n — Eo)(n\ = H — Eq, the commutation relation (H — £'o)p!f r o = i(VV") Po, 
an integration by parts, and Gauss' Law (S7 2 V — 47rZ^ i S^ 3 '(ri)), the denominator V is easily evaluated: 

V = ( * | p ■ (H - E Q )p | W ) = 2nZ W 2 (0), 

but the logarithmic factor makes the numerator (3 much harder to evaluate. Even for a very simple one-electron 
system such as the hydrogen atom, (3 cannot be evaluated in closed form, though several rapidly convergent methods 
can be used to evaluate it to high accuracy |l0| , |^T| , and the Bethe logarithm of the electronic ground state 
of was recently evaluated numerically ]l3|. For a two-electron system such as the helium atom, whose unknown 
wavefunction 'J'o must be represented by an expansion in a large basis set, the numerical challenges are even more 
daunting. 

In the early 1960's C. Schwartz recast the numerator as integral over the virtual photon energy k 

/? = lim ( -K (W | p p | W )+Vln(K)+ 

K—>oo \ 
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and thereby replaced the insuperable difficulties associated with accurately summing over an infinite number of bound 
and continuum eigenstates of H with the more tractable difficulty of numerically integrating an accurate representation 
of the matrix element of the resolvent (H — Eq + k)~ x for small, intermediate and large values of k. When k is very 
large, Schwartz found it sufficient to approximate the matrix element with a simple asymptotic formula. For smaller 
values of k, the action of the resolvent is solved explicitly as the solution of a system of linear equations in a suitable 
basis with p-wave symmetry. For intermediate k the convergence was greatly improved by including a single function 
which has the same leading-order asymptotic behavior as the true solution as k — > oo. 

Despite growing problems with the numerical linear dependence of his basis as the number of basis functions was 
increased, Schwartz was able to compute for the 1 1 S ground state of the neutral helium atom a Bethe logarithm 
of 4.370(4) Rydbergs, which yielded a theoretical ionization potential for this state in agreement with the best 
experimental values available at that time, and which remained unsurpassed until very recently. 

The results presented in this letter were generated by an approach very similar to that used by Schwartz, in which 
the integral in Eq. ([!]) is split into a low k region /3l and a high k region /3h- The counterterms in Eq. (Q) are then 
brought inside the integral to cancel explicitly the divergent behavior at large k: 
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where rph(k) and ifin(k) are solutions of the equations 

(H - E Q + k) Mk) = P*b, (3) 
(H - E Q + k) Vta(fc) = (H - E )p* . (4) 

Since H possesses overall rotational symmetry, the solutions iph(k) and ip-a(k) have a total angular momentum 
quantum number which can differ by only ±1 from that of Wo- In this work <Po has ^-symmetry, so i/^h(k) and ipn(k) 
have P-symmetry. 

An elegant derivation of Eq. (Q) can be found in the work of Forrey and Hill jl2| , which examines Schwartz's method 
from a fresh perspective and provides many useful computational techniques. We evaluate the two integrals in Eq. 
(0) numerically, using the procedure described by Forrey and Hill, computing the matrix element of the resolvent at 



each integration knot by solving variationally for -0l or tpn in Eq. (||) and Eq. (||). When k is very large, we use the 
asymptotic approximation 14 
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The constants C and D have been computed in closed form only for the hydrogen atom; in this work they are estimated 
by extrapolating the values generated by the solution of Eq. (Q) at successive integration knots. This equation was 
solved explicitly at each successive knot, running in the direction of increasing k, until the relative difference between 
successive extrapolated estimates of C was roughly 1%. For larger k the resulting asymptotic formula was used. For 
the helium ground state our estimates of C and D are 4.988(1) and -18.8(3) respectively, with the errors resulting 
mainly from extrapolation uncertainty. These estimates can be compared with the value 5.18 computed by Schwartz 
Q for C and the value -20±3 he assumed for D. 

The non-relativistic wavefunction •Z'o was computed variationally using our modification fuj, II of the basis 
set first developed by Frankowski and Pekeris which exploits knowledge of the analytic structure of the true 
wavefunction at the 2- and 3-particlc coalescences to improve the convergence of the variational trial function to the 
exact unknown wavefunction: 

% = ^ c "(^( s > *i u ) ± 4>u{s, -t, u)) 
<j)u{s,t,u) = s n t l u m {\nsye- as+ct 

where s, t, and u are the Hylleraas coordinates defined by s = ri+r2, t = r2—r\ and u = r\i and the ± sign is chosen 
so that the product of <Po and the spin function is antisymmetric under exchange of the electrons. 

Our bases for representing iph{k) and ipn(k) include functions of four different types. The fc-independent functions 



X\P = n(j> v (s,t,u) ± r 2 (j>u{s,-t,u) 

together with the single function x^ = P&o provide a good solution space for small k. 

For large k the solution ipn(k) becomes concentrated in fc-dependent regions of configuration space for which one 
electron is very close to the nucleus and the other electron is much further away, so it is essential to use explicitly 
fc-dependent basis functions. Of primary importance is the 'Schwartz function' x^ 3 \ an approximate solution of Eq. 
that reproduces the first two terms in the asymptotic expansion in Eq. (Q): 

To help approximate that part of ipn(k) which is orthogonal to the 'Schwartz function', we also use a fourth set of 
functions x^\ which are symmetrized sums of products of single- variable Laguerre functions Ct(Rj) = Li(Rj)e~ Rj I 2 
of the three perimetric coordinates i?i =ri+r 2 — ri 2 , R2—Ti^f2 J rfi2 1 and R% = — ri+r2 + ri2: 

xl 4) = (ri C p {aR x ) C q {bR 2 ) C r (cR 3 )) ± (n ~ r 2 ). 

Combinations of the exponential parameters a, b, and c can be chosen to reflect the strong 'in-out' correlation in tpji(k) 
for large k. For any k the overlap matrix elements for these basis functions are very small or zero far from the main 
diagonal, which enables us to avoid the severe problems with numerical linear dependence which prevented Schwartz 
from using a large basis of functions of the form of powers of r 2 , r 12 times a highly asymmetrical exponential 
of r\ and r 2 . We set a = b + c to eliminate from x^ an y exponential ri 2 -dependence, which would complicate the 
evaluation of matrix elements between these functions and the other types of basis functions. Analytic considerations 
jUI suggest that the integrand is optimised if b w (2k) 1 ' ' 4 and c m Z. We coarsely search the parameter space in the 
neighborhood of these values of b and c seeking to maximize the two integrands of Eq. (^|) in accordance with the 
variational principle described in fl^| . 

The calculation of /?l was fast and straightforward. In this case k was small enough that there was no need to 
include explicit fc-dependence in the basis, x^ was omitted altogether, and a single average value of the parameter 
b was used in the x^ functions, independent of the value of k at a particular integration knot. We solved for ipL{k) 



in a basis with 92 Xv functions, the x^ function, and 120 Xu functions. The parameters b and c were varied to 
maximize the integrand. Changes in the integrand due to small variations in b and c were used to assess convergence. 

The /?h integral was computationally expensive, primarily because including the 'Schwartz function' x^ requires 
evaluating algorithmically complicated matrix elements. Since x is intended primarily to accelerate the convergence 
for very large k, and since over half of the knots in our integration scheme correspond to k < 40, we chose to omit 
-^(3) f rom the basis for knots below k w 40. At each node we solved for ipn in a basis consisting of 92 X" functions, 
the x^ 2 ' function, the x^ function (for high k), and 220 x^ functions. We then recomputed the solution of Eq. (^) 
after first reducing the number N of x^ functions in the existing matrices to study convergence of the integrand. A 
simple polynomial fit in the variable 1/N was applied to the sequence of results with N — 220, 165, 120, and 84 to 
generate the values of the Bethe logarithms in this letter. The error associated with the finiteness of the basis for ipn 
was taken as the entire difference between the extrapolated value and the value corresponding to N — 220. 

Other sources of numerical error arise from the numerical integration itself (for which there are good analytic error 
bounds p2|), and the finiteness of the basis used to approximate \Po. The latter error is assumed to be comparable 
to the relative error in T> in all cases. For neutral helium, independent runs with less accurate representations of 
<?b indicate that this estimate of this error is somewhat conservative. The numerical integration was parametrized 
to keep the absolute error in /3h and /3l below 10 -8 . The results of independent calculations carried out for neutral 
helium with a coarser mesh were consistent with the analytic error bound. 

The uncertainties assigned to the Bethe logarithms in this letter are the sums of the uncertainty due to extrapolation 
of ipn and the uncertainty due to approximation of $o in a finite basis. The uncertainty in tpi, and the numerical 
integration error bounds are negligible by comparison. 

Our Bethe logarithms for the 1 1 S, 2 1 S, and 2 3 S states are listed in Tables |, ||, and [II, respectively. The values of 



ko have been divided by Z 2 to illustrate their approach to the hydrogenic limit as Z becomes large. Scaled values of 
the nonrelativistic binding energy E U r/Z 2 and V/Z 4 are also listed to provide some measure of the accuracy of 'Pq. 
Uncertainties in T> were computed by comparison with highly accurate results for <8{r\) > provided by Drake Jl6| |. 

The exact hydrogenic limits of ln(fco/Ry) pill , p7f| and of -Ehr and V are displayed in the bottom row of each table, 
labeled by l/oo (exact). Immediately above the bottom row, in the row labeled l/oo, we list the hydrogenic values 
and the corresponding uncertainties computed using the method described in this letter with \/Z = so that the 
I/V12 term is removed from the Hamiltonian. 

For Z = 1 the Hamiltonian H has a single bound state of 1 1 S symmetry. As Z — ► 1 from above all the singly 
excited bound states of a two-electron ion disappear into the continuum as the 'outer' electron moves infinitely far 
away. Hence as Z — > 1 from above, the energies and all other finite-range properties of the states should tend toward 
those for a single hydrogen atom in its ground state with Z = 1. The approach of the Bethe logarithms and other 
properties toward their hydrogenic values as Z — > 1 from above is visible in Tables II and III for the 2 1 S and 2 3 S 
states, respectively. 

We have fit our ionic results to the 1 jZ expansion developed by Goldman and Drake [H 

ln(fc /Ry) = Co + d /Z + C 2 /Z 2 + ■ ■ ■ 

C = ln2 + 21nZ + ln(fc ff /Ry), (6) 

where ln(fcjj/Ry) is the weighted sum of the two hydrogenic Bethe logarithms corresponding to the state. Table [i"v| 
displays the results of a three parameter polynomial fit for C\, C%, and C3 using data for Z = 4, 8, and 16. The listed 
uncertainties come from the formal propagation of error through the regression formula and do not include truncation 
errors from higher-order terms in the expansion. 

Our results for the 1 1 S state and the 2 1 S states of neutral helium are in complete agreement with the recent 
calculations of Korobov and Korobov p9| . The most accurate previous value of the Bethe logarithm of the 2 3 S state 
came from Goldman and Drake's lst-order 1/Z expansion ]l7[ , pq . A numerical comparison of results for neutral 
helium appears in Table 0. 

Preliminary values of our Bethe logarithms for the 1 1 S, 2 1 S I pO l and 2 3 S states of helium were used in a recent 
comparison of theory and experiment by Drake and Martin |21f|. The values in this letter make slight corrections 
to the theoretical ionization energies of the 2 1 S and the 2 3 S levels in that work, while the 1 1 S state is unaffected. 
Modifying "Bethe log cor." contribution in Drake and Martin's Table II to include the values in this letter yields the 
theoretical results in Table fVl| , which are compared with results from several recent experiments [^2| , |2^] . |24j , |^5) , 
§§• 
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TABLE I. ln(fc /(Z 2 Ry)), Enr/Z 2 and V/Z A for the 1 X S state. 



l/Z 


ln(fc /(Z 2 Ry)) 


Enr/Z 2 (a.u.) 


V/Z* (a.u.) 


1/1 


2.992 97(5) 


-0.527 751015 308 


2.067 80(4) 


1/2 


2.983 864(2) 


-0.725 931094 259 


2.843 815 67(5) 


1/3 


2.982 624(2) 


-0.808 879 268 074 


3.189 069 9(1) 


1/4 


2.982 503(1) 


-0.853 472 889 901 


3.376 853 2(1) 


1/8 


2.982 948(2) 


-0.924 321798 793 


3.677 270 31(4) 


1/16 


2.983 448(1) 


-0.961551275 290 


3.835 858 39(2) 


l/oo 


2.984128 6(7) 


-1.0 


4.0 


l/oo (exact) 


2.984128 556 


-1.0 


4.0 



TABLE II. ln(fc /(Z 2 Ry)) , Enr/Z 2 and V/Z* for the 2 X S state. 



l/Z ln(fc /(Z 2 Ry)) Enr/Z 2 (a.u.) V/Z 4 (a.u.) 

1/1 (exact limit) 2.984128556 -0.5 2.0 

1/2 2.980115(1) -0.536 493 511514 2.056 896 21(2) 

1/3 2.976362(2) -0.560 097416177 2.103163 60(3) 

1/4 2.973 976(1) -0.574 054 618 459 2.132 593 60(7) 

1/8 2.969 797(4) -0.597 793 082 931 2.185 583(1) 

1/16 2.967459(3) -0.610 956 015 708 2.2163203(6) 

l/oo 2.964 977 7(4) -0.625 2.25 

l/oo (exact) 2.964 977 593 -0.625 2.25 



TABLE III. ln(fco/(Z 2 Ry)) , Enr/Z 2 and V/Z* for the 2 3 S state. 



l/Z ln(fc /(Z 2 Ry)) E n r/Z 2 (a.u.) V/Z A (a.u.) 

1/1 (exact limit) 2.984128 556 -0.5 2.0 

1/2 2.977 742(1) -0.543 807 344 559 2.074 008 93(2) 

1/3 2.973852(1) -0.567 858 596 952 2.124 087184(9) 

1/4 2.971735(1) -0.581072 911861 2.152 566 566(2) 

1/8 2.968414(2) -0.602 260114 376 2.1991479(4) 

1/16 2.966 705(1) -0.613 440 895 056 2.224 062 6(2) 

l/oo 2.964 9776(2) -0.625 2.25 

l/oo (exact) 2.964 977 593 -0.625 2.25 



TABLE IV. Coefficients of the l/Z expansion using a 3 parameter fit. The exact values of Ci are due to Drake p7j . 

l/Z coeff 1_^S 2^S 2^S 

Ci (exact) -0.0123 03(1) 0.040 771(1) 0.027 760(1) 

Ci -0.0123 2(5) 0.040 78(10) 0.027 73(5) 

C 2 0.0228(8) -0.016(2) -0.0010(6) 

C 3 0.002(2) -0.011(6) -0.007(2) 



TABLE V. Comparison of ln(fco/Ry) for neutral helium. The uncertainty in the lst-order l/Z expansion due to uncalculated 
higher-order terms could not be readily estimated until more exact calculations were done. 



State 


lst-order l/Z expansion 


Schwartz 


Korobov 


This work 


1 X S 


4.364(?) 


4.370(4) 


4.370 157 9(5) 


4.370 159(2) 


2 X S 


4.372(?) 




4.366 4091(5) 


4.366 409(1) 


2 3 S 


4.365(?) 






4.364 036(1) 



TABLE VI. Ionization potentials, in MHz, as described in |2lJ, but with theoretical values corrected slightly by the results 
in this letter, a: from | p2| , b: from p3| , c: an average of values from |24j and [^BJ, d: from pq |. 



State This work Experiment Difference 



1 L S 


5 945 204 226(91) 


5 945 204 238(45) a 
5 945 204 356(48) 6 


12(102) 






130(103) 


2 l S 


960 332 040.9(25.0) 


960 332 041.01(15) c 
1 152 842 742.87(6) d 


0.1(25.0) 


2 3 S 


1 152 842 738.2(25.2) 


4.7(25.2) 



